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Abstract 

We  have  considered  some  peculiarities  of  light  reflection  by  several  types  of  fractal 
surfaces.  The  two-dimensional  Weierstrass  function  for  modeling  of  fractal  surfaces  was 
used.  The  computer  simulation  of  real  rough  surfaces  was  performed  by  means  of  this 
function.  The  electromagnetic  waves  scattering  indicatrix  was  obtained  for  the  concrete 
fractal  surfaces  based  on  the  scalar  Kirchhoff  theory.  The  computer  simulation  of  light 
reflection  by  fractal  surfaces  was  performed. 


1.  Introduction 

All  the  surfaces  are  rough  in  some  way  or  another.  Therefore  the  studying  of  the  scattering  of 
electromagnetic  waves  by  rough  surfaces  is  an  important  and  interesting  theoretical  and  experimental 
task.  The  investigation  of  the  peculiarities  of  such  scattering  is  an  important  first  of  all  for  a non- 
destructive control  of  real  surfaces.  Real  surfaces  are  most  adequately  described  by  fractal  functions 
as  they  are  neither  pure  deterministic  nor  pure  casual.  Last  time  the  fractal  nature  of  a number  of 
different  surfaces  have  been  experimentally  determined  (sea  surface  and  sea  bottom,  Earth  relief, 
cloud  surfaces,  thin  film  surfaces,  deposited  on  the  substrate,  etc.). 

The  aim  of  this  work  was  the  simulation  of  the  electromagnetic  waves  scattering  indicatrix  based 
on  the  scalar  Kirchhoff  theory  for  particular  surfaces.  The  analogues  calculations  were  performed  by 
others  authors  [2],  but  our  results  have  some  distinctive  features,  particularly  in  the  expression  for  the 
averaged  scattering  coefficient  there  are  some  additional  terms  which  can  significantly  influence  on 
the  resulting  scattering  indicatrix  at  a certain  geometry  of  the  experiment. 


2.  Theory 

We  have  chosen  the  two-dimensional  Weierstrass  function  z(x,  y)  for  the  simulation  of  rough 
surfaces 
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were  cw  is  a normalizing  factor;  q > 1 is  the  fundamental  spatial  frequency;  K is  the  fundamental 
wavenumber  of  the  surface;  D is  the  fractal  dimension  ( 2 < D < 3 );  N,M  are  numbers  of  tones;  tynm 
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is  a phase  term.  The  example  of  a rough  surface  simulated  using  the  function  (1)  is  shown  on  Fig. 
1(a). 

Let  us  consider  the  wave  falling  on  the  rough  surface  S at  an  angle  ©j  and  scattering  at  a polar 
angle  ©2  and  at  an  azimuth  angle  03.  We  will  be  finding  the  scattered  field Es(r,t)  based  on  the 
scalar  Kirchhoff  method  [4]: 

Es  (r ) = -ikrF(Ql , 02©3  J exp  [ikty(x0  , y0  )]dx0dy0  + Ee(r ) , (2) 

*0 

where  k is  the  wavenumber  of  the  incident  wave, 


F(0j,©2,03)  = +B2  +C2) 


is  the  angle  factor,  R is  the  reflection  coefficient, 

A = sin  0j  - sin  02  cos  ©3 , 
B = -sin02sin03, 


C = —cos  0[  - cos  02 , 
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is  the  phase  function, 
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is  the  bound  term, 
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The  above  formalism  is  valid  under  the  following  conditions [4]:  the  incident  wave  is  monochromatic 
and  plane;  the  scattering  surface  is  rough  inside  a certain  square  and  smooth  outside  its;  the  surface 
dimensions  are  much  large  than  the  incident  wavelength;  all  the  surface  points  have  finite  gradient; 
the  reflection  coefficient  is  a constant  across  the  surface  area;  the  scattered  field  is  observed  far  from 
the  surface. 

After  some  transformations  from  Eq.(2)  considering  Eq.(l)  we  obtain  the  expression  for  the 
average  scattering  coefficient 


h 


(4) 
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where  (ls)  = (EsE*},  I0  = 


(2kXYcos&t  Y 


Kr 


is  the  intensity  of  a wave  reflected  from  the  respective 


smooth  surface.  Neglecting  the  terms  higher  than  = kcwCq^D  3^“),  this  expression  has  the 

following  approximate  form 
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where,  cis  the  root  mean  square  height  of  the  surface  roughness,  X and  Y are  the  dimensions 
reflecting  area 
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3.  Numerical  results 

We  have  calculated  the  averege  reflected  coefficient  as  a function  of  02  and  03  (the  scattering 
indicatrix)  based  on  Eq.  (5).  We  have  assumed  that  R - 1 , in  other  words  we  have  not  considered  the 
real  dependence  of  the  reflectance  R on  the  wavelength  X and  on  the  incident  angle  ©j.  The 
example  of  the  scattering  indicatrix  is  shown  on  Fig.  1 (b). 


4.  Conclusion 

The  analysis  of  the  obtained  results  leads  to  some  inferences: 

• The  waves  scattering  is  symmetrical  relatively  the  plane  of  incidence; 

• The  most  intensity  of  the  scattered  waves  is  observed  in  the  specular  direction; 

• There  are  others  directions,  where  some  splashes  of  intensity  are  observed; 

• The  picture  of  the  reflection  complicates  with  increasing  of  the  surface  large  scale 
homogeneity. 

These  peculiarities  are  due  to  combination  of  chaotic  character  and  self-similarity  of  the  real  surface 
relief. 
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Fig.  1 a)  The  simulation  of  the  fractal  surface  by  means  of  Weierstrass  function  K = 6.3,  N = M = 5, 
D = 2.5;  b)  The  reflection  coefficient  log(p^)  for  the  fractal  surface  with  D = 2.5,  q = 1.8; 
N = M=  10. 
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